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A METHOD FOR NUMERICALLY CALCULATING THE AREA AND 
DISTRIBUTION OF WATER IMPINGEMENT ON THE 
LEADING EDGE OF AN AIRFOIL IN A CLOUD 
By Norman R. Bergrun 


SUMMARY • ■ 

A method is presented for determining, "by step— by— step integra- 
tion, the trajectories of water drops around any body in two- 
dimensional flow for which the streamline velocity components are 
known or can be computed. The method is general and considers the. 
deviation of the water drops from. Stokes 5 law because of speed and 
drop size. 

The equations are presented In general- form and then, to 
illustrate the procedure, water-drop trajectories are calculated 
about a 12— percent— thick symmetrical Joukowski profile chosen to 
simulate an NACA 0012 section. 

The method provides a means for the relatively rapid calculation 
of the trajectory of a single drop without the utilization of a 
differential analyzer. 

In addition, consideration Is given to the maximum possible rate . 
of water-drop impingement on a body. 


INTRODUCTION 

As part of a comprehensive research program directed toward a 
fundamental understanding of the mechanism of thermal ice— prevention 
for airplanes, the NACA has undertaken an investigation of the 
performance of a heated wing or empennage section in specified icing 
conditions. One of the first essentials of such an investigation is 
a method for estimating or calculating the area over which water will 
strike the wing, and the distribution of water impingement over that 
area. 
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Previous Investigations of this nature have considered the 
water-drop trajectories about cylinders (references - ! and 2) and 
those about an airfoil (reference 2). In both of these studies the 
assumption was made that the water drops were sufficiently small to 
obey Stokes' law of resistance. A more recent investigation for 
cylinders (reference 3) shows that, for the velocities of airplanes 
and the drop sizes frequently present in clouds, Stokes' law does 
not apply and the force acting on the drop can be determined only, 
from a knowledge of the drag coefficient for spheres. 

In reference 3, therefore, the methods of previous investigators 
are modified to include more accurate values of the drop drag coeffi- 
cient. Equations are established for the determination of the drop 
trajectories around cylinders, spheres, and ribbons. These equations 
are derived in forms suitable for use in a differential analyzer and, 
by utilizing an analyzer, the trajectories are calculated and plotted. 

The designer of a heated wing, desiring to know the rate and 
area of water Impingement on the leading edge in a specified cloud at 
a given flight speed, might assume that the impingement will bo 
identical to that for a cylinder with a radius equal to the wing 
leading-edge radius. There is some question, however, as to the 
accuracy of this assumption for the larger drop sizes and for wing 
sections with small leading-edge radii. The designer Is therefore 
confronted with the desirability of employing a computation method, 
preferably without the necessity of a differential analyzer, which 
will provide, some Indication of the area and distribution of water 
Impingement. 

In this report, the water-drop trajectory equations used in 
reference 3 have been modified to establish a step-by— step Integra— 
tion method applicable to any two-dimensional flow for which tho 
streamline velocity components are known or can be approximated. If 
drop trajectories for a large range of airspeeds and drop diameters 
are required, the computation time Is large, and - the desirability of 
access to a differential analyzer becomes evident. The integration 
method presented herein, however, does permit tho calculation of 
any desired number of trajectories without resort to, an analyzer; and 
it also provides a means for estimating the error which will be 
incurred by replacing the airfoil by a cylinder with a radius equal 
to the leading-edge radius of the airfoil.' In addition, water- 
Impingement-’data over the entire airfoil surface can be obtained by 
tho integration method. 
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SYMBOLS 


Coefficients 


radius of a drop, microns (3*28 x lCT e ft) 

average acceleration of a drop over an interval of time, 
feet per second per second 


acceleration of a drop at the "beginning of an interval 
of time, feet per second per second 


dimensionless ratio 



chord of airfoil, feet 
dimensionless concentration: factor 
drag coefficient for spheres 

efficiency of deposition, of water drops on a body 

dimensionless drop^-inertia quantity, defined in equation (k) 

dimensionless drop— inertia quantity for which no drops 
impinge upon a body 


L unit length of span of an airfoil 

M intensity of water-drop impingement at a point on the surface 

of a body, pounds per hour, square foot 

Mp rate of water-drop impingement on a body, pounds per hour, 

foot span 

p nondimensional distance 

R local Reynolds number of a drop relative to the air stream 

% Reynolds number of a drop moving with free— stream velocity 

in still air — 

s distance along the surface of a body from zero percent 

chord, feet 


dimensionless streamline parameter 
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t unit of time^jj ^ , seconds 

u x , u y* u 2 velocity components of air, relative to free— stream 
velocity 


U free— stream velocity of air, feet per second 

v Xj Vy, Vz velocity components of a water drop, relative to free— 
stream velocity 


v 

w 

x, y, z 

ymax 

a 

(3 

A 

e 

9,t 

Xq 

P 

p 

Pi 

* 

Subscripts 


Average velocity of a drop over an interval, of time, 
feet per second 

liquid water content of the air, grams per cubic meter 
(62.4xicr 6 ,ib/cu ft) 

rectangular coordinates, feet 

maximum ordinate of airfoil 

condition for which no drops will strike a body 

proportionality factor 

increment 

thickness parameter for transforming a oircle into a 
Joukowski airfoil 

coordinates of a circle of which the origin is at the 
center 

range of a drop in still air, feet 

absolute viscosity of air, pounds mass per second;, foot 

1 

density of a water drop, pounds per cubic foot 
density of air, pounds per cubic foot 
stream function 


denotes starting conditions where time is counted as zero 


o 
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1 

denotes 

cr 

denote s 

n 

denotes 

n+i 

denotes 


on© interval of time 
a critical condition 
nth interval of time 
interval of time of (n+i) 


ANALYSIS 

To obtain area and concentration of water— drop impingement at 
the leading edge of an airfoil, the path of the drops as affected 
by the air stream must be traced. In order to compute these paths, 
or trajectories, it is necessary to define the motion of a drop 
through the air. The equations which depict this motion, and which 
have been modified to include the effect of a drop departing from. 
Stokes* law of resistance have, been shown by Langmuir and Blodgett in 
reference 3 ho take the following form: 


Kv x ^ „ M 
dx 2 k 


(Tt-ux). 


(1) 


y dy 


Odr 

nr 


(vy-Uy) 


(2) 


dx 

dt 


= — Vx 


dy _ 

~ 


(^Ru^f ~ ^ v x ^x^ 2 + ^y^y^ 2 


(3) 


The dimensionless quantity K, used by Glauert (reference 2 ), 
is evaluated by the equation 


k = - ^ 

9 \Pi J 


/& \ g / Piffc 
\c ) r p 


) 


- w 
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where in a system of consistent units , p is density of the drop. 
Pi is density of air, a is radius of the drop, IT is free—stream 
•velocity, (i is absolute viscosity of air, and c is the chord of 
the airfoil. 


It can be seen from equations (l) and (2) that 

1 


dv x 

"dT 


CpH 

24K 


(vx— Ux) 


dv y 

dt 


CdR 

24K 


(Vy-Uy) 


Since the average velocity of a drop over a small interval of time 
may be expressed by 


v = Y n+i+vn 
2 


and the average acceleration by 


a = a n+i +a n 
2 


the velocity and position of a drop at the end of any interval of 
time may be expressed by the following equations: 


(n+x) Vx n 24K ^ 




(5) 


CdR 


r y (n+i) r ?n 24K u yn ^ 


( 6 ) 



NACA TN Ho. 1397 


7 


x ( n+i ) = Vxn At — i (v Xn — u Xn ) At 2 + Xn (7) 

2 |_ 2kK 


y(n+x) “ v 7n At — — |jj| (v yn -Uy n ) ] At 2 + j n (8) 

where the subscript n denotes a point at the beginning of some 
interval of time, and the subscript n+i denotes a point at the 
ena of that interval of time. The main assumption made in these 
equations is that the acceleration of the drop is constant over the 
interval of time selected. 

Equations ( 3 ), ( 5 ), ( 6 ), ( 7 ), and ( 8 ) may "be used to calculate 
the trajectory of a drop by a step-by-step process, provided that 
a sufficient range in values of air-velocity components u x , u y , and 
the drag coefficient function CpE/24 are available. When Stokes f 
law of resistance holds, the value of CpH/2h is unity. The 
quantities u x , Uy, and CpR/Sk, are best used in large-scale plots 
to insure a fair degree of accuracy in the computations. In addition 
to the^large-scale plots required, certain assumptions must be made 
in regard to the conditions of a drop at the start of its trajectory. 
This information is needed in order to continue with the step-by-step 
computations. It appears that the simplest manner in which to start 
the computations is by assuming the velocity of the drop to be the 
same as that of the air at a fairly large distance forward of the 
airfoil, and then estimating the velocity differentials Vx-u x and 
vy—uy between the air and the drop at a new point An estimation of 
the magnitude of these differences may be made by the following 
equations which are derived in appendix A: 



/ 
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If y 0 is the ordinate of a trajectory for a given drop size at 
infinity, the quantity of water per unit time and area which progresses 
toward a body is 


UwLy 0 

Since this same amount of water distributes itself, over the 
leading edge of the body, then, 

UvLy 0 = MLb 


and for a small increment 

M = Uw ( 11 ) 

ds 

Since dy 0 /ds is proportional to the concentration of water impinge- 
ment at various points on a body in a water cloud, it can be appro- 
priately called a concentration factor. Thus, 

M. =s UwC ' ' “ ' (12) 


The total water collected by a body is the summation of all the 
degrees of concentration of water-drop impingement over the leading 
edge of a body for a unit length of span such that 


n 8 

Mp = / M&s = TJwy ^y E (13) 

J o 

In equation ( 13 ), E .is efficiency of water disposition over the 
leading edge, and is defined in terms of the maximum ordinate of the 
body, so that 


y °limit 

^max 


(HO 


where yon^t. is the oxdinate of a starting trajectory which Just 
impinges upon the airfoil, and y^a* is the maximum ordinate of 
the body. 


When drops are quite small, they have insufficient momentum to 
carry them to the body, hence, there is a critical size of drop which 



FACA ITT No. 1397 


9 


possesses just sufficient energy to impinge upon the body. This 
critical condition is given by 


a = ^-cr 


du x _ 1 
dx _ ¥ 


M15) 


To determine the critical drop size (i.e., the critical K value) , 
it is only necessary to evaluate the quantity du^/dx at the 
stagnation point of the airfoil. Because the quantity du x /dx is 
a property of the streamlines around a body at the stagnation point, 
it has dif ferent values dependent upon the body chosen. 


CCMPOTATIOF PKOCEDUKE 

Hie example given herein to explain the computation procedure 
is for a symmetrical Joukowski airfoil, 12 percent thick, at 0° angle 
of attack. This airfoil was chosen because the streamline functions 
could be calculated with relative ease, and because its profile 
closely simulated that of an KA.CA 0012 Airfoil. Trajectory data 
for an FACA 0012 section were needed. in connection with thermal ice- 
prevention heat— transfer computations . - - 

Streamlines were obtained about the Joukowski airfoil for a 
region estimated to cover the paths of drops of particular interest. 

The equations used to define the contour of the streamlines and the 
values of velocity components along streamlines were derived from. 
Glauert (reference 4) . These equations are given in appendix C. 

Streamlines- calculated .for the Joukowski airfoil are shown in 
figure 1. In. addition, the contour of the NAGA 0012 airfoil is 
presented to show it3 comparison with the Joukowski profile over the 
region where it was expected that trajectories would hit for drop 
sizes typically found in nature. 

Air—' velocity— component data were used in the forms shown in 
figures 2(a), 2(b), 2(c), 2(d), and 3. Figure 2(a) gives the 
x— component of air velocity in terms of the x— coordinate in the 
vicinity of. the airfoil, while figure 2(b) presents the same parameters 
several chord lengths forward of the airfoil. Figure 2(c) depicts the 
y— component of velocity in terms of the y— coordinate in a region near 
the airfoil. However, because of their small magnitudes, y-components 


x Derived in appendix B. 
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of velocity well forward of the airfoil had to he treated somewhat 
differently than x— velocity components in order to obtain values . 
sufficiently accurate.* It was found that the y— velocity components 
showed a good degree of change with respect to a change in ratios of 
the coordinates at any point. It was therefore found desirable to 
plot the y— velocity components in the manner shown in figures 2 (d) and 
3 . Both x— and y— velocity components were plotted so that their values 
could bo obtained to four significant figures. 

Because a water-drop trajeotory can cross streamlines,, and 
because air-velocity components are conveniently plotted In terms of 
streamlines, it is imperative to be able to interpolate between the 
streamlines to ascertain on what streamline a drop is located at a 
given time. Interpolation between streamlines can be accomplished by 
plotting points on streamlines in terms of ratios of the point 
coordinates. • Figure 4 shows the relationship obtained when this was 
done. The simplicity of this method of interpolation lay in the fact 
that the streamlines displayed a logarithmic separation, which permits 
the interpolation to be carried on by laying a logarithmic scale 
properly across the streamlines. It should be observed that, as a 
trajectory approaches the ordinate axis (taken In the c Deputations 
to originate at the stagnation point), the ratio of the coordinates 
asymptotically approaches infinity. Therefore, in performing the 
step— by— step computations, it is advisable to abandon this Interpolar- 
tion procedure when values of the coordinate ratio are of the order 
of 1200 percent. For small values of x/c forward of the airfoil, the 
value of the stream function can be obtained for the drop with the aid 
of the information given in figure 1 . 

Values of the drop drag-coefficient function Cj)R/2k and R 
were obtained from reference 1 , and plotted to a suitable scale for 
use in the computations. These values are given in table I. As the 
stejv-by—step process progresses, new values of R and CpR/24 are 
selected from the plots. 

From equations ( 3 ), (5), ( 6 ), ( 7 ), and ( 8 ) it is apparent that 
trajectories can be calculated for both ranges of K and R-q values, 
where Ru is defined by ■ 


% = 


2 aUpi 


However, because of limited time available for performing the 
calculations, a set of trajectories (corresponding to several K 
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values)' was computed for only an Eg- value of 95*65- The value \6f 
Eg- was ’ chosen "because it corresponded to certain airspeed, altitude, * 
and dropr-size conditions of particular interest; ' - - 

Starting coordinates of trajectories -were arbitrarily chosen for 
selected K values.’ As an example, by means of equations ( 9 ) and (10), 
for a K value of O.O 58 I and starting coordinates Xq= 210 percent 
chord, yo = 0*750 percent chord, it was determined that 

O.OOOl 


v yi"“ u 3 r i - 


- .= 0.9979 


v yi = 1.5X1CT 5 - 


CpB _ 

IT ” 


1.001 


Eg - 95*65 


A sample calculation for obtaining these values is. shown in 
appendix D. With these values as starting values, a continuation 
of the computations is shown in table II for the trajectory. This 
table expresses equations ( 3 ), (5)j (6), (7)* and (8) in a convenient 
computational form and is used in conjunction with figures 1, 2(a), 
2(b), 2(c), 2(d), 3, and k. 


EESULTS 

From calculations to determine the trajectories of water drops, . 
the area and distribution of water impingement over the_ leading 
edgs of an airfoil can be obtained. In addition, the total amount of 
water per unit length of span cm Ido calculated. 
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If trajectories for one Ey value are plotted in the manner 
shown in figure 5, the distance along the surface to the stagnation 
point included by a trajectory may he measured. Then by talcing 
included distances from a large-scale plot similar to figure & 
figure like that shown in figure 6 is. obtained. The limiting trajec- 
tories beyond which no drops, will impinge upon the airfoil surface 
may be obtained .by connecting the end points for the several computed 
trajectories. This is shown by a broken line in figure 6 . It is 
interesting to note in figure, 6 , that for ’the case of the 

problem reduces to one of plotting the airfoil ordinates against the 
included surface distances. 

By obtaining the slopes of the curves in figure 6 , according to ll 
equation (ll), the concentration factor C may be obtained at jj 

different points over the leading edge of the airfoil. This operation \[ 
was performed graphically, and the results are shown in figure 7 . \ 

The ratio of the limit starting ordinates of trajectories 
y°limit ' fco " tlle max i ;nium ordinate of -the airfoil defines the efficiency 
of drop impingement (equation (it)). Therefore, when such ratios are 
obtained by means of figure 6 , and plotted with respect to their 
corresponding K values, an indication of the efficiency of water 
deposition on an airfoil can be obtained. The value of K corre- 
sponding to zero impingement was estimated to bo approximately 
0.0009 by means of equation (15)- The value of du x /dx needed 
to estimate this value of K was determined graphically by using a 
large plot of figure 2(a) . A presentation of the efficiency of 
water-drop impingement on the Joukowski airfoil, for a value of 
Ru=95*65, * s Bilo ' wn in figure 8 . 

As a matter of Interest, the paths of several drops of different 
inertia but the same starting points were obtained in relation to 
the velocity field about the airfoil. These relations are shown in 
figures 9 (a.) and 9 (b)., and were obtained by noting from the computa- 
tions (as performed by the procedure suggested in table H) the 
cor re sponding air— velocity components at different times during a 
drop trajectory. 


DISCUSSION 

The procedure demonstrated herein for calculating water-drop 
trajectories about a body in two-dimensional flow is convenient to 
use when aid from a differential analyzer is not available. In the 
use of the method, however, certain errors are' involved. One of the 
errors Is Involved in the ability to predict, with reasonable accuracy. 
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the starting values of the velocity differentials v Xl — Ux x and 
Vyi~ u yi for small drop sizes (K values less than approximately 
0.02). It appears that, in order to calculate trajectories for these 
small drops, either the air— velocity components must he obtained to • 
great accuracy or else a better method of approximating 'the ^starting ~ 
values of the differentials' must be derived. Another error lies in 
the accumulation, during successive integration steps, of errors 
involved by the assumption of an average constant acceleration over 
the time interval selected. Near the airfoil, errors accumulate more 
rapidly than they do well forward of the airfoil because of the 
increased rate at .which values of the air— velocity components change. 
(See figs 9{a.) and 9(b).) Large changes in airr-velocity components 
mean large changes in selected values of 0]jR/2b which multiply the . 
velocity differentials vx-ux and vy— Uy. Therefore, it is desirable 
to reduce the magnitude of the arbitrary time interval .when a 
trajectory approaches the leading edge. To determine a suitable time 
interval to apply in the computation of a trajectory as it approaches 
a body, a criterion is to permit the drop to travel, in the time 
interval selected, a distance in which the air— velocity x— component 
values do not change by more than approximately 10 percent. 

To greatly reduce the amount of work involved in computing a 
set of trajectories, starting ordinates of trajectories can be taken 
very close to the body, and time intervals can be chosen as large 
as is commensurate with the criterion for their selection. When this 
was done for the computation given in table II (the s tar ting ordinate 
being taken at half a chord length forward of the airfoil), area of 
water-drop impingement was reduced 15 percent. To start the calcula- 
tions so near the airfoil, it was assumed the drop followed the 
streamline up to the starting point. Because of this assumption the 
starting ordinate was large by approximately 10 percent, compared to 
2 percent for the computations given in table II. , 

With the set of trajectories calculated by the step-by— step 
process and presented in figure 6, a comparison can be made between 
the impingement of water drops on the Joukowski airfoil and its 
equivalent cylinder for a sms.ll rang9 in drop sizes. Such a comparison 
is made in figure 10, where the airfoil was used as the basis for 
comparison. For area of impingement, the equivalent cylinder can be 
used without much error to a drop^-size diameter of _ approximately 20 
microns. However, for- rate of impingement, appreciable difference 
is not encountered until drop diameters are above 15 microns/ It Is 
of interest to note that meteorological data obtained In flight 
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(reference 5 ) indicate a mean-effective drop— size diameter 2 for 
cumulus clouds of 15 to 25 microns. Furthermore, mean-effective drop 
diameters as high as 50 microns have "been measured in clouds of the 
altostratus category. Since, in accordance with the definition of 
mean— effective drop diameter, half the water contained in a cloud . , 
sample is composed of drops larger than the H©an-effective diameter, 
it is apparent that the use of an equivalent cylinder-may "be subject 
to considerable error when applied to typical icing conditions in 
which large drop-size diameters have been observed. 

The method for computing trajectories of water drop's in two- 
dimensional flow can be extended to the case of three-dimensional 
■flow about a body of any shape . . To extend the method, the trajectory 
of a drop must be expressed for the third plane. The additional 
trajectory equations, needed are: 


Kv, 


dv z PdE 

Tz = ~ mr (vz ~ Uz) 


( 16 ) 



In addition, equation ( 3 ) which connects the variable drag coefficient 
Cp with the local Reynolds number of the drop at a given time becomes 


(^2^) = (V;£-U x ). 2 + (vy— Uy } 2 + (v z — u z ) 2 (17) 


Thus 'when equation (l 6) is placed in forms similar to equations ( 5 ), 
(6), .(7), and (8) in the analysis, and equation (17) is solved in a 

2 Mean-effective drop diameter as used in the present report is the 
sizo of drop in a cloud sample for which the amount of liquid water, 
existing in water drops larger than that drop is equal to the amount 
of liquid water existing in drops smaller than the drop. 
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maimer similar to that indicated in table II, a three-dimensional case 
may he handl ed by Btep— by— step integration when air— velocity components 
and starting conditions are known. Application of trajectory computa*- 
tions for three— d ime nsional flow can be made to the case of windshield 
ice prevention, where it is desirable to know the rate a£ which water 
impinges upon a windshield in determining heat requirements to keep 
a windshield clear of ice. Equations (l6) and (17) are basic ones 
which can be adapted also -for use with a differential analyzer. 

Ihe concentration factor C in equation (12) for the three- 
dimensional case is now evaluated by 



(18) 


In equation (l8), A c is an area which is outlined by several- 
starting trajectories; and Aq is the area, on the surface of the 
body, which is outlined when these trajectories intercept the body. 

. The rate of water impingement Mp in equation (12) is evaluated 
for the three-dimensional case by 


Mp = UwEAfrontal (19) 


In equation (19), E, the efficiency of water-drop deposition on the 
body, is now expressed by 


_ Aon mi t 
Afrontal 


whore is the area included by starting trajectories which 

Just impinge upon the body, and Af r ontal is the frontal area of the 
body. 
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As' has already been mentioned, the limit distribution of water- 
drop impingement on the leading edge of a body c&n be obtained from 
the body ordinates. It is interesting to observe that, for very high 
speeds, the limit condition .K=» -is approached. Hence, it is 
believed that a comparison -of the water-drop impingement over several 
different airfoils for the limit condition is of interest. Such a 
comparison is made in figure 11. For these limit conditions, . any body 
having a leading— edge radius will have a value of the concentration 
factor at. the leading edge of unity; whereas, any body without a 
leading-edge radius will' have a much lower value of the concentration 
factor depending upon the slope of the contour near the stagnation 
point. According to equation ( 13 ), the area under the curves in 
figure 11 is proportional to rate of icing per unit span, and in 
accord with theory, for airfoils with the same value of maximum 
thickness, regardless of location, the areas under the curves in 
figure 11 are equal. The difference between oach of the throe 
airfoils lies in the manner in which. water drops distribute thomsolvos 
over the airfoil. The reason for this difference is because a limit- 
ing water-impingement distribution corresponds to the case whoro drops 
are either so large or are traveling so fast that they are undeflectod 
by forces in the air stream. Therefore, the effect of shifting the 
location of maximum thickness Is to change the degree and position of 
water-drop impingement. Thus the limit water-drop-impingement distri- 
bution for a given body can be made to have any desired profile by the 
proper combination of contour and location of maximum thicknoss. 

CONCLUSIONS 

From the method outlined in this report to obtain trajectories 
of water drops forward of a body, the following conclusions can be 
drawn: 


1. By a step— by— step process which considers deviations from 
Stokes* law of re si stance, the area and distribution of water-drop 
impingement on a body of any shape and angle of attack can be calculated 
when the air— velocity components are known. 

2. The rate and area of water-drop impingement on an NACA 0012 

airfoil having an 8— foot chord Is approximated with reasonable 
accuracy by the impingement on a cylinder of radius equal to the 
leading— edge radius of the airfoil for only drop diameters below 20 
microns at spoods of about 200 miles per hour and at an altitude of 
7000 feet. - 
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3. The critical drop-size diameter corresponding to zero impinge 
ment can “be determined for a two-dimensional tody if the velocity 
components are known at the stagnation point. 

b. The area and distribution of water-drop impingement on an 
airfoil in motion through a cloud of known drop-size distribution c an 
bo obtained by calculating the area and distribution of impingement 
for each drop— size range present. 

5* Thin airfoils of different sections (such as wedge— shape, 

NACA 0006, and circular— arc airfoils), having the same exposed frontal 
areas, will have approximately the same rates of water-drop impingement 
at high speeds, but the distribution, of water-drop impingement will 
be different for each section, and will depend upon the rate of change 
of the slope of the airfoil surface and the location of maximum 
thickness. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. ( July 1947. 
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APPENDIX A 

Derivation of Starting-Point Equations 

At a position Xq, y 0 well forward of the airfoil, it is 
assumed that a drop is moving with a free— stream velocity u^. 


U X 0 - T Xo 
u yo = v 7o 


The position of the drop after an interval, of time At is 


xi = x 0 + Ax 
yi » y Q + Ay 


where Ay may he expressed in terms of Ax hy 


Ay 


u yo 

u x 0 


(Ax) 


The velocity of a drop, starting initially from a point 
after an interval of time At is 


vxi = + Av x 



+ Av y 


Then 

(Al) 


(A2) 
o* 7o-> 


(A3) 
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From 'basic trajectory equations 


vx dv x = B (yx~^x) dx 


(Ah) 


where 


B = 


CpR 

2kZ 


By substituting equations (Al) and (A 3 ) into equation (A^) 


Vxi ( 7 n-u x ) - B (v Xl -u Xl ) Ax 


and by solving for v XlJ 


v Xx = 


(u Xq +BAx) ± a / (u 3 Co +BAx) a -ABu 2;i Ax 


(A5) 


In a manner similar tp that by which equation (A5) was derived, 
the initial value of the y— component of the drop velocity v v may be 
obtained. The value of vy^ becomes 


- 


- (BAy-Uy 0 ) (BAy-Uy Q )2+4-Bu yi Ay 


(A 6) 


2 
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APPENDIX B 

Derivation of Equations for the Condition Where 
No Drops Impinge Upon a Body 

In deriving the relationships which describe the limiting value 
of K below which no drops will impinge upon a body, it is necessary 
to consider the basic trajectory equations and the flow field at the 
stagnation point of the body. 

Near the stagnation point' of a body, the velocity of both air and 
of drops which nearly follow the streamlines is so small that CjjR/2k 
may be taken as unity. In addition, the air-velocity components are 
nearly linear in relation to their distance from the stagnation point. 
Glauert shows in reference 2 that K is proportional to the range of 
a drop in still air; and in reference 3> Langmuir denotes this range 
by Xq, where 



c 


Since K is a function of the air— velocity components and 
because, near the stagnation point of the body x-ccmponents of air 
velocity; according to reference 3? may "be simplif ied to 


ux » ap (Bi) 


and the range of a drop rqay be written in terms of a distance 
p forward of the stagnation point of a body. This distance is 
defined by 


P = ~ (B2) 

A. s 

Equation (l) may be written 


Kv* 


avx 

dx 




vx-ap 


(B3> 
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To solve equation (B3), v x may be replaced in manner similar 
to Ug. by 


v x = £p 


(B4) 


where by differentiating 


dv £ s pdp + pd (3 


Then 


df = P+P dt (®5) 

By differentiating equation (B2) with respeot to x, 

Ep (B6) 

and hy substituting equations (B4) and (B 6 ) into equation (B 3 ), 




(B7) 


By setting equations (B5) and (B7) equal to each other and rearranging 

*C033Zl&j 


= __ ■ {3dg 
5 «H3+(3 £ 


(BS) 


gives 


For the integration of equation (B 8 ), Piero e’s tables (reference 6 ) 


/: 


BdJ3 ZnX 


_ _g Tdg 

OVi J V 


co-p+P*' 2h 2h u X 

4 ' 


(B9) 
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where 


X = a+gP+hp 2 


ana 



f 2h.P+gA 

V-/7 J 


In equation (B9) 


e = ikxh-g 2 

Since g=>~l and h=l, equation (B8) becomes 


- Znp =-i-ln(ctr-P+p 2 ) - ,, 1 t an-i ( 

J 4(i-l ■ \ 


In equation (BIO) , it is apparent that 
in evaluating the function of p. 



Is a critical 


A substitution of equation (Bl) into equation (B2) gives 


. Ku x = cue 


and by differentiating' with respect to x;' 


du.,. 

Kcr -5 — = ct, 
dx. 


(BIC) 

value 


Since equation (310) is dependent 'solely upon the trajectory 
equations, the value of a may always be taken as l/b in evaluating 
the critical K value below which no drops will impinge on a body. 
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However, tho quantity dux/dx is dependent upon the shape of a "body., 
the attitude of the 1 ody in the air stream, and the type of flow around 
the body. It, therefore, must be evaluated at the stagnation point for 
the flow arcuni the body under consideration. 


• APPENDIX C 

Equations for Obtaining Streamlines and Telocity 
Components about a Symmetrical Joukowski 
Airfoil , 12 Percent Thick 

The contour of a streamline is given by 
x = (r cos 3 + £) ! 1 


1 '\ 

r z+eZ+ 2 re cos 9 J 


. . / 1 

y = (r sin €>)( 1 

\ r 2 +e a +2re cos B 


) 


in which 


r 


s 


2 sin Q 


+ /( ■-- -W l-fg) 1 
V \2 sin & 


(Cl) 


In equation (Cl) 


s = ~ 


t 

U 


where ^ , a stream function, can have any desired value. In addition, 
for a symmetrical Joukowski airfoil, 12 percent thick, the value of e 
is 0.1021. 
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The components of air velocity in terms of free— stream velocity 
are given "by ■ 


U x = U3.U2 - V1V2 


Ujr = UlV2 + uavi 


where 


, - (^ £ ) 

ui 1- ~ — cos 2 0 


r 

,2 


(1+e) 

Vl » sin 2 6 


u 2 = 


" pCp - 1 )^ 2 1 

. (:d-i) 2 +f 2 j 


(C2) 


vp. 


7 


p -1 

(d-i) 2 +f 2 j 


( 03 ) 


In equations (C2) and (C3)j 


D = T z cos 20 + 2 re cos 0 v + £‘ 


F = 2r sin 0 (r cos 0 + ’£ ) 
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APPENDIX D 

Computation of Starting Values for the 
Step-tv— Step Computations. 

According to equation (A5) , the x— velocity component of a drop is 


v Zl 


(u Xq ^x) 2 -43u Xi Az 


In order to obtain starting values with this equation, this equation 
must be solved by a trial— and— error method . A procedure for a trial— 
and-error solution Is demonstrated for the starting values used in 
table II- 

Por the conditions of table II, ' 


Xo = 2.1: y D = 0.0075 
K = 0 . 0581 ; Ey = 95-65 


Several chord lengths forward of the airfoil, as a first approx- 
imation, the quantity vy x — Uj x may be considered zero. Thus 


R~±BU (v Xl -u Xl ) 


Assuming a value of— (v Xl — u Xl ) = 0.0001, the value of R is 


E ~ 0.00956. 
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for which the value of CjjR/ 24 is 


CdR 

a 1.001 


With x Q as 2.1 chord lengths forward of the airfoil., and x* as 
2.0, the corresponding magnitudes of air-velocity components are 
u x =0-9979 and u Xl =0.9978- The value of Ax is then 0.1, and the 
value of B is 17-228, from which the value of v Xl from equation 
(A5) is 0-9979- Therefore, the initial assumption of 
— (yxi'”U 3 Ci) =0.0001 was very nearly correct. 

Since the quantity (vy x ~ u yi ) contribute a very little to the 

initial value of Cjfi/2k 3 once the approximate value of v Xl — Uxi 
has been determined, an approximation of v y;L — 1 u yi may be made by 
equation (a 6) without changing the initial value chosen for CdR/ 24- 
For the trajectory calculated in table II, vy x — uy x was assumed to be 
0.1x10“ 5 , and it was also found by equation (A2) and (A6) to be of 
this magnitude. 

After the differentials v x — and v y 1 ~ u y JL are evaluated, 

they are entered in the appropriate columns in table II, and the 
process continued by the steps indicated in the table. 
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TABLE I 
CdR 

VALUES OF AS FUNCTIONS OF R 


R 

CdR 

24 

R 

cdr 

24 

0,00 

1.00 

200 

6.52 

•05 

1.009 

250 

7.38 

•1 

1.013 

300 

8.26 

•2 

1.037 

350 

9.00 

•4 

1.073 

400 

9.82 

•6 

1.103 

500 

11.46 

•8 

1.142 

600 

12, 97 

1.0 

1.176 

800 

15.81 

1.2 

1.201 

1000 

18.62 

1.4 

1.225 

1200 

21.3 

1.6 

1.248 

1400 

24.0 

1.8 

1.267 

1600 

26.9 

2.0 

1.285 

1800 

29.8 

2.5 

1.332 

2000 

32.7 

3.0 

1.374 

2500 

40.4 

3.5 

1.412 

3000 

47.8 

4. 

1.447 

3500 

55.6 

5. 

1.513 

4000 

63.7 

6. 

1.572 

5000 

80.0 

8. 

1.678 

6000 

96.8 

10. 

1.782 

8000 

130.6 

12. 

1.901 

10000 

166.3 

14. 

2.009 

12000 

204. 

16. 

2.109 

14000 

243. 

18. 

2.198 

16000 

285. 

20. 

2.291 

18000 

325. 

25. 

2.489 

20000 

365. 

30. 

2.673 

25000 

470. 

35. 

2.851 

30000 

574. 

40. 

3.013 

35000 

674. 

50. 

3.327 

40000 

778. 

60. 

3.60 

60000 

980. 

80. 

4.11 

60000 

1175. 

100. 

4.59 

80000 

1552. 

120. 

5.01 

10? 

1905. 

140. 

5.40 

1.2X10* 

2234. 

160. 

5.76 

1.4X10? 

2549. 

180. 

6.16 

1.6X10 

2851. 
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IABLB II.- COMPUTATION OF A WATER-DROP TRAJECTORY 
[ RU ® = 9.15 x 10 3 i K = 0.0681] 


3 

a 

K 

► 

d 

X 

d 

i|s 

o |n 

# 

d 

i 

► 

^ d 
H 

? 

c 

H 

0.06 

0.9979 

0.9978 

1.001 

17.2289 

0.0001 

0.001722 

0.997814 

.06 

,9978 

.9976 

1.002 

17.2461 

.0002 

.002828 

.997673 

.10 

.9977 

.9976 

1.002 

17.2461 

.0002 

.002983 

.997376 

.10 

.9974 

.9972 

1.003 

17.2633 

.0002 

.005884 

.996987 

.10 

.9970 

.9968 

1.003 

17.2633 

.0002 

.003228 

.996664 

.10 

.9967 

.9964 

1.004 

17.2805 

.0003 

.004562 

.996208 

.10 

.9962 

.9960 

1.003 

17.2633 

.0002 

.003691 

.995849 

.10 

.9968 

.9966 

1.006 

17.3160 

.0003 

.006042 

.996246 

.10 

.9962 

.9960 

1.004 

17.2805 

.0002 

.004233 

.994822 

.10 

.9948 

.9944 

1.008 

17.3494 

.0004 

.008188 

.994003 

.10 

.9940 

.9935 

1.009 

17.3666 

.0005 

.008735 

.993130 

.10 

.9931 

.9922 

1.016 

17.4871 

.0009 

.016263 

.991504 

.10 

.9916 

.9910 

1.009 

17.3666 

.0006 

.008762 

.990629 

-.10 

.9906 

.9892 

1.024 

17.6248 

.0014 

.025168 

.988112 

.10 

.9881 

.9870 

1.019 

17.6387 

.0011 

.019503 

.986162 

.10 

.9861 

.9840 

1.038 

17.8667 

.0022 

.038625 

.982299 

.10 

.9823 

.9796 

1.060 

18.0723 

.0028 

.060684 

.977241 

.10 

.9772 

.9725 

1.083 

18.6403 

.0047 

.088373 

.968404 

.10 

.9684 

.9610 

1.126 

19.5804 

.0074 

.143 512 

.954063 

.10 

.9640 

.9400 

1.220 

20.9983 

.0141 

.295089 

. 924544 

.01 

.9245 

.8925 

1.379 

23.7349 

.0320 

.760561 

.916938 

.01 

.9169 

.8845 

1.382 

23.7866 

.0324 

.771690 

, .909222 

.01 

.9092 

.8750 

1.396 

24.0276 

.0342 

.822269 

.901199 

.01 

.9011 

.8650 

1.410 

24.2685 

.0362 

.878495 

.892414 

.01 

.8924 

♦8530 

1.432 

24.6472 

.0394 

.971445 

.882700 

.01 

.8827 

.8380 

1.468 

25.2668 

.0447 

1.12943 

.871406 

.01 

.8714 

.8220 

1.498 

25.7831 

.0494 

1.27384 

.868668 

.01 

.8587 

.8020 

1.544 

26.5749 

.0667 

1.50595 

.843609 

.01 

.8436 

.7790 

1.688 

27.3322 

.0647 

1.76591 

.825950 

.01 

.8260 

.7510 

1.648 

28.3649 

.0760 

1.97605 

.806190 

.01 

.8062 

.7140 

1.738 

29.9139 

.0922 

2.76776 

.778612 

.01 

.7786 

.6680 

1.870 

32.1859 

.1106 

3.56015 

.743011 

.01 

.7430 

.6110 

1.990 

34.2513 

.1320 

4.52165 

.697796 

.01 

.6978 

.642 0 

2.180 

37.5215 

.1658 

6.84570 

.639339 

.01 

.6393 

.4720 

2.360 

40.4476 

.1673 

6.76844 

.571656 

.01 

.6717 

.4400 

2.600 

44.7504 

.1317 

5.89161 

.612739 

.01 

.6127 

.5440 

2.536 

43.6489 

-.0413 

1.80100 

.530749 

.01 

.5307 

.7000 

2.640 

43.7177 

-.1693 

7.39926 

.604742 
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TAELB II— COHTIEPSD. CCMFUTATIC* 0? A EATER-DROP TRAJJCTORT 
[ E^a - 9.16 x 10 *i K - 0.05811 


> 

« 

$ 

1 

I 

aw 

O JCM 

a 

t 

i 

N 

y 

R 

1 * 

H 

x/o 

m 

+ 

o 

X 

/ 

? 

o 

H 

X 

/ 

0.049895 

0.0025 

0.00086 

2.1600 

0.049893 

1.96011 

1.6000 

1.60 

.049890 

.0026 

.00141 

3.6350 

.049887 

1.90022 

1.6861 

1.76 

.0997673 

.0100 

.00160 

14.9179 

.099752 

1.80047 

1.7274 

1.90 

.0997376 

.0100 

.00194 

19.4212 

.099718 

1.70076 

2.0250 

2.20 

.0996987 

.0100 

.00161 

16.1412 

.099674 

1.60108 

2.3270 

2.50 

.0996664 

.0100 

.00228 

22.8102 

.099644 

1.60144 

2.6256 

2.90 

.0996208 

.0100 

.00180 

17.9669 

.099603 

1.401 B 4 

3.0997 

3.50 

.0995849 

.0100 

.00302 

30,2147 

.099565 

1.50229 

3.7941 

4.20 

.0995245 

.0100 

.00211 

21.1686 

.099605 

1.20279 

4.4969 

5.0 

.0994822 

.0100 

.00409 

40.9446 

.099441 

1.10336 

5.3662 

6.0 

.0994003 

.0100 

.00436 

43.6770 

.099357 

1.00399 

6.4667 

7.4 

.0993130 

.0100 

.00813 

81.3160 

.099232 

.90476 

8.0882 

9.2 

. 0991604 

.0100 

.00437 

43.7638 

.099107 

.80665 

10.032 

12.0 

.0990629 

.0100 

.01258 

126.841 

.098937 

.70671 

13.452 

16.6 

.0988112 

.0100 

.00976 

97.615 

.098714 

.60800 

17.061 

20.0 

.0986162 

.0100 

.01931 

195.128 

.098423 

.50958 

22.216 

29.0 

.0982299 

.0100 

.02629 

252.922 

.097977 

.41160 

34.555 

43.0 

.0977241 

.0100 

.04419 

441.868 

.097282 

.31432 

49.996 

70.0 

.0968404 

.0100 

.07176 

717.660 

.096123 

.21820 

87.284 

120.0 

.0964053 

.0100 

.14754 

1476.44 

.093930 

.12427 

160.689 

256.0 

.0092454 

.0001 

.38026 

58.0280 

.0092074 

.11506 

369.724 

660.0 

.0091693 

.0001 

.38580 

38.6795 

.0091308 

.10693 

438.621 

760.0 

.0090922 

.0001 

.41113 

41.1134 

.0090531 

.09688 

612.687 

880.0 

.0090119 

.0001 

.43924 

43.9247 

.0089680 

.08791 

600.934 

1020.0 

.0089241 

.0001 

.48572 

48.5722 

.0088766 

.07903 

702.636 

1200.0 

.0088270 

.0001 

.66471 

56.4710 

.0087705 

.07026 

825.222 

1600.0 

.0087140 

.0001 

.63692 

63.6920 

.0086603 

.06161 

996.717 

1800.0 

.0086866 

.0001 

.76297 

76.2970 

.00851158 

.06310 

1203.09 

2300.0 

.0084560 

.0001 

.88295 

88.2950 

.0083478 

.04476 

1464.59 

2800.0 

.0082695 

.0001 

.98802 

98. e 020 

.0081607 

.03669 

1862.32 

3800.0 

.0080619 

.0001 

1.37888 

137.888 

.00802301 

.02867 

2411.94 

4800.0 

.0077861 

.0001 

1.78007 

178.007 

.00760812 

.02096 

3126.30 

7000.0 

.0074301 

.0001 

2.26077 

226.077 

.00720403 

.01376 

4378.08 

10000.0 

.0069779 

.0001 

2.92286 

292.285 

.00668568 

.00707 

6300.37 

16000.0 

.0063933 

.0001 

2.38422 

338.422 

.00606497 

.00102 

9939.82 

26000.0 

.0067165 

.0001 

2.94580 

294.580 

.00542197 

-.00440 

16031.3 

42300.0 

.0061273 

.0001 

.90050 

90.060 

.00603754 

-.00944 

27786.8 

66000.0 

.0063074 

.0001 

3.69963 

369.963 

.00493763 

-.01438 

59665.1 

61700.0 
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.93425 

1.11178 

1.9676 

2.0472 

2.9390 

6.7844 

8.6666 

20.0041 

32.7159 

104.311 

290.276 

311.379 

367.313 

419.066 

497.366 

674.778 

804.283 

1096.910 

1345.410 

1937.68 

2388.06 

3873.70 

6626.92 

9699.63 

15060.2 
26269.0 

27213.2 
22034.9 


1.72289 
2.82594 
2.97599 
3.02038 
2.98617 
4.74177 
6.94330 
7.02799 
8.69366 
10.9947 
16.2247 
19.4418 
34.2044 
36.0816 
51.6462 
121.208 
156.626 
372.882 
654.047 
2190.36 
6889.67 
7406.66 
8826.61 
101.701>&£ 
122.687X10 
170.495 XLtf 
207.369XKL 
291. 603 XH? 
367.730>l CT 
549.621XLC? 
714.362x10? 
124.678*10 
192.729>4Cf 
363.946Xlcr 
609.1 47X1CT 
1175.548X30 
1187.826)40? 
963.316X10 


1.68614 
1.72744 
2.02504 
2.32708 
2.62560 
3.09978 
3.79411 
4.49691 
5.56626 
6.46676 
8.08622 
10.0324 
13.4628 
17.0610 
22.2166 
34.3334 
49.9969 
87.2841 
160.689 
369.724 
458.621 
612.687 
600. 934 
702.636 
826.222 
996.717 
1203.09 
1494.69 
1862.32 
2411.94 
3126.30 
4373.08 
6300.37 
9939,82 

16031.3 
27786.8 
39665.1 

49298.3 



4.999 

8.728 

16.068 

3.697 

4.386 

6.126 

6.009 

7.026 

6.262 

9.957 

12.030 

14.545 

16.623 

24.119 

31.263 

43.730 

63.003 

99.398 

160.313 

277.668 

396.661 


8 

6 

19 

22 

26 

29 

34 

41 

49 

69 

73 

91 

117 

153 

196 

283 

422 

681 

1190 

2602 

404 

476 

557 

652 

764 

910 

1099 

1349 

1638 

2137 

2769 

3750 

5337 

8120 

12986 

21909 

33726 

44482 


76.608 

76.616 

76.636 

76.667 

76.682 

76.711 

75.745 

76.786 

75.836 

76.894 

75.967 

76.068 

76.176 

76.328 

76.624 

76.807 

77.229 

77.910 

79.100 

81.702 

82.106 

82.582 

83.139 

83.791 

84.565 

85.465 

86.664 

87.913 

89.561 

91.688 

94.457 

98.207 

103.544 

111.664 

124.650 

146.659 

180.285 


.003877 

.003979 

.004200 

.004448 

.004727 

.006043 

.006403 

.006819 

.006306 

.006878 

.007666 

008406 

.009455 

.01080 

.01269 

.01607 

.01876 

•02479 

.05626 

.066716 

.07136 



96.66 

96.65 

95.66 

96.66 

96.65 

95.66 

96.66 

96.65 

95.65 

96.65 



.9110 

1.6941 

12.450 

3.3320 

1.9100 
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TABLE II.- COUC LOT ED. COHFUTATIO* OF A RATER-DROP TSAJECTCBT 
[ Ed* . 9*16 X 10*| I ■ 0.0681] 


















ORD/NJATE N£ , PERCENT 


^ JOUKOWSKI PROFILE 
- NACA 0012 PROFILE 


-STREAMLINE 
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distance: along chord lime Y c , PERCENT 
STREAMLINE FIELD ABOUT A SYMMETRICAL JOUKOWSKi AIRFOIL 12 
PERCENT THICK SHOWING A COMPARISON BETWEEN THE SURFACE 
STREAMLINES FOR AN N/CA 0012 AND THE JOUKOWSKI AIRFOIL 





X-COMPONEUT OF AIR VELOCITY U K1 PERCENT OF FREE STREAM 



DISTANCE ALONG CHORD DATUM LINE, >£ , PERCENT 


FIGURE 2A~ AIR- VELOCITY -COMPONENT PROFILES FOR VARIOUS STREAMLINES ABOUT A SYMMETRICAL JOUKOWSKl 
AIRFOIL |2 PERCENT THICK.. X-VEUDCfTY COMPONENTS IN PROXIMITY TO THE AIRFOIL 



NACA TN No. 1397 Fig. 2a 



FIGURE 2B.- AIR-VELOCITY- COMPONENT PROFILES TOR VARIOUS STREAMLINES 
ABOUT A SYMMETRICAL JOUKOWSKI AIRFOIL \l PERCENT THICK. 
X-VELOCITY COMPONENTS SEVERN. CHORD LENGTHS FORWARD OF THE AIRFOIL 


X-COMPlOMENT OF A\R VELOCITY, U*, PERCENT OF FREE5TREAM 



A68I 'ON Ni VOVN 


qS '3M 




NACA TN No. 1397 


Fig. 2c 
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DISTANCE FROM CHORD DATUM UKJE.,%, PERCENT 
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FIGURE 2Cr AIR-VELOCITY-COMPONENT PROFILES FOR VARIOUS 
STREAMLINES ABOUT A SYMMETRICAL JOUKOWSKl 
AIRFOIL 12 PERCENT THICK. Y- VELOCITY COMPONENTS 
IN PROXIMITY TO THE AIRFOIL 




Y-COMPOH&NT OF All? VEIjOCITY, U Y , PERCENT OF FREE STREAM 



I Z M 6 a \0 20 HO £0 40 Wo 20 0 w 

RATIO OF Y- COORDINATE TO Vf- COORDINATE, IN PERCENT 


FIGURE 2D.' AIR- VELOCITY COMPONENT PROFILES FOR VARIOUS STREAMLINES ABOUT A SYMMETRICAL 
JOUKOWSK1 AIRFOIL 12 PERCENT THICK. Y- VEUXlTY COMPONENTS IN TERMS OF 
coordinate RATIOS 


. 2d NACA TN No. 1397 



STREAMLINE VALUE, 5 

FIGURE 3.- COORDINATE RATIOS OF POINTS ON STREAMLINES AS A FUNCTION OF Y” COMPONENTS 
OF AIR VELOCITY FOR. A SYMMETRICAL JOUKOWSKI AIRFOIL 12 PERCENT THICK. 


RATIO OF Y-COORDINATE TO X-COORDINATE - \N PERCENT 



8 *3T£ 


A68T '°N NX VOVN 



RATIO OF Y- COORD! NATE TO X- COORDINATE IN PERCENT 



.1 Z H & & 1.0 2.0 40 LX) SjO lao 20.0 


DISTANCE ALONG CHORD DATUM UNL IN PERCENT 

FIGURE 4r STREAMLINES ABOUT k SYMMETRICAL JOUKCWSKI AIRFOIL 12 PERCENT THICK 
AS A FUNCTION OF THEIR COORDINATES. 


Fig. 4 NACA TN No. 1397 




RATIO OF Y-COORD1NATE TO X-COORDfNATE,^, IN PERCENT 


to 
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FIGURE 4 .-C 0 NCLUDED 


NACA TN No. 1397 Fig. 4 cone. 





FIGURE 5.- RELATION OF WATER-DROP TRAJECTORIES TO STREAMLINE 
FIELD ABOUT A SYMMETRICAL JOUKOWSK! AIRFOIL 12 
PERCENT THICK. 


Fig. 5 NACA TN No. 1397 
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o 4 8 12 16 20 24 28 i 

DISTANCE ALONG SURFACE OF AlftFaL FROM LEAPING EDGE, % IN PERCENT 
FIGURE 6 - AREA OF DROP INTERCEPTION CORRESPONDING TO DROPS OF DIFFERENT INERTIA 
FOR A SYMMETRICAL JOUKDWSKI AIRFOIL |Z PERCENT THICK. 




FREE STREAM REYNOLD’S 
NUMBER OF DROP, Ru* 95.65 


\ DIMFNaONLfSf 



a^^DROP-lMEKTlA QUANTITY. K- 
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6 8 

DISTANCE ALONG SURFACE OF AIRFOIL FROM LEADING EDGE.,% I hi PERCENT 


FIGURE 7 ~ DISTRIBUTION OF WATER-PROP IMPINGEMENT OVER A SYMMETRICAL JOUKOWSK\ AIRFOIL 
12. PERCENT THIOL- 
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EFFICIENCY OF DROP IMPINGEMENT, E, IN PERCENT 

figure a - water- drop -collection efficiency tor a sti-imetrical joikowski airfoil 

12 PERCENT THICK., 00 


NACA TN No, 1397 Fig. 



X-COMPONENT OF AIR VELOCITY, U K , PERCENT OF FREE STREAM 



~6 -4 -2 o a 4 6 8 10 

DISTANCE ALONG CHORD DATUM LINE, >Sfc, PERCENT 


FIGURE %•- INFLUENCE OF VELOCITY FIELD CM THE PATH OF WATER DROPS HAVING DIFFERENT INERTIA 
FOR A SYMMETRICAL J0UKOW9KI AIRTOlL 12 PERCENT THICK. X-VELDCnY COMPONENTS 


Fig. 9a NACA TN No. 1397 



Y-COMP0NENT OF A'R VELOCITY, U Yt PERCENT OF FREE STRBAM 


NACA TN No. 1397 


Fig. 9b 
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DISTANCE FROM CHORD DATUM LINE, PERCENT 
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FIGURE 9B- INFLUENCE OF VELOCITY FIELD ON THE. PATH OF WATER 
DROPS HAVING. DIFFERENT INERTIA FOR A SYMMETRICAL 
J0UKDW5K1 AIRFOIL 12 PERCENT THICK. Y- VELOCITY 
COMPONENTS 



TUAN ON THE EQUIVALENT CYLINDER 



DROP-SR& DIAMETER, MICRONS 

FIGURE 1 0 .- COMPARISON OF RATE. AND AREA OF U/AT&RrDROP IMPINGEMENT ON A 

SYMMETRICAL JOUICOWSVCI AlfcFon- \2 PERCENT THICK AND A CYLINDER. 
WITH RADIUS EQUAL TO THE AIRFOIL LEADING EDGE RADIUS. 



!A TN No. 1397 
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FIGURE II.- DISTRIBUTION OF WATER. DROPS ON THREE DIFFERENT AIRFOILS WHEN THE. 
DROPS ARE UN DEFLECTED BY FORCES IN THE AIR. STREAK. 


NACA TN No. 1397 


